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ABSTRACT 


A  finite  element  computer  code,  STAGS  C-1  and 

holographic  Interferometry  were  used  to  determine  the 

effects  of  Interior  cutouts  (2x2,  2x4,  and  4x4  Inch) 

on  the  first  five  natural  frequencies  and  mode  shapes  of 

curved  Graphite  Epoxy  panels.  The  panels  are  a  quasl- 

Isotroplc  layup  [0,  -45,  45,  90]  with  a  12  inch  chord  and 

s 

height.  Both  the  finite  element  and  holographic  analysis 
were  conducted  using  c 1 amped-c 1 amped  boundary  conditions. 

The  vibration  branch  of  STAGS  C-1  1 s  an  energy 
technique  based  on  small  displacements  and  linear  elastic 
stress-strain  relationships.  When  compared  with  the  time 
averaged  holograms  of  the  experimentally  determined  natural 
frequencies  and  mode  shapes,  the  two  techniques  show  a  close 
correlation  of  both  frequency  and  shape. 

It  was  found  that  for  the  2x2  Inch  cutout,  the  mode 
shapes  change  very  little  while  the  natural  frequencies 
displayed  a  small  decrease  for  the  higher  modes.  The  2x4 
Inch  cutout  retained  the  general  mode  shape  of  the  solid 
panel  for  the  first  two  modes.  The  third  through  the  fifth 


modes  shapes  were  changed  by  this  cutout  and  the  loss  of 
panel  stiffness  was  visible  .  The  4x4  inch  cutout  exhibit 
both  a  switch  in  symmetry  of  the  first  two  modes  and  a 
general  decrease  in  natural  frequencies. 


I.  INTRODUCTION 


Background 

Composite  materials  are  used  frequently  as  both  primary 
and  secondary  structural  elements  in  the  construction  of 
modern  aircraft.  This  trend  will  probably  accelerate  as  the 
higher  performance  levels  demanded  in  military  aircraft  and 
the  increased  economy  sought  by  civilian  aviation  can  be 
achieved  through  extensive  use  of  these  high  strength, 
lightweight  materials.  In  order  to  use  composite  materials 
more  extensively  in  the  aircraft  industry,  a  complete 
understanding  of  dynamic  response  is  indicated.  Curved 
shells  with  rectangular  cutouts  were  chosen  because  this 
shape  Is  often  found  in  aerospace  structures  and  no  closed 
form  solution  currently  exists.  Various  studies  of  a 
simillar  type  have  been  conducted  on  flat  plates.  Ashton 
and  Whitney  [9]  were  able  to  form  a  solution  for  symmetric 
orthotropic  plates  where  no  coupling  existed  between 
extension  and  bending.  Simlliariy,  Raju  [10]  was  able  to 
derive  exact  solutions  for  circular  isotropic  plates  with 
central  circular  holes.  Jones  and  Morgan  [11]  were  able  to 
derive  exact  numerical  solutions  for  buckling  and  vibration 
problems  for  snt 1 symmetr i c  cross  ply  shells.  Their  results 


, 
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show  that  the  effect  of  coupling  between  bending  and 
extension  on  buckling  loads  and  vibration  frequencies  dies 


out  rapidly  as  the  number  of  layers  Increases.  However,  for 
unsymmetrlc  shells  the  coupling  dies  out  very  slowly.  And 
In  a  si ml  liar  study,  Monahan  [12]  applied  a  finite  element 
solution  to  Isotropic  rectangular  plates  with  arbitrary 
cutouts  In  order  to  determine  their  natural  frequencies  and 
mode  shapes.  Results  were  compared  with  the  results  of  an 
ho  I ograph 1 c  ana  lysis. 

The  logical  extension  of  these  efforts  was  to  extend 
the  analysis  to  a  more  complex  structure  and  shape.  By 
choosing  curved  panels  with  their  frequent  use  In  the 
aerospace  field.  It  Is  hoped  that  a  valuable  contribution 
can  be  made  to  the  technology  base  of  composite  materials. 
Additionally,  a  preliminary  evaluation  of  the  STAGS  Cl 
Finite  Element  Code  In  modelling  real  panels  Is  desired. 


APPROACH 


A  parallel  effort  to  utilizing  the  Finite  Element 
Method  (FEM)  and  laser  holography  was  undertaken  to 
determine  the  first  five  natural  frequencies  and  mode  shapes 
of  graphite  epoxy  shells  with  various  retangular  cutouts. 


The  STAGS-Cl  computer  code  developed  Oy  B.  0.  Almroth,  F.  A. 
Brogan,  and  G.  M.  Stanley  of  the  Lockheed  Palo  Alto  Research 
Laboratory  was  chosen  due  tc  its  ease  of  application  to 
orthotropic  shells  and  user  friendliness.  Highly  accurate 
time  averaged  holography  was  used  to  experimentally 
determine  the  natural  frequencies  and  mode  shapes.  An 
outstanding  characteristic  Inherent  to  this  technique  Is 
that  Its  data  output  Is  a  photographic  Image  readily 
compared  with  the  contour  plots  generated  from  FEIi 
calculations.  The  natural  frequencies  and  modes  shapes 
extracted  by  the  two  methods  were  compared  In  order  to 
analyze  the  effectiveness  of  the  FEM  to  predict  the 
experimental  results. 


V 


The  governing  equation  of  motion  for  a  body 
experiencing  free  vibration  Is  derived  from  Newton's  Second 


The  Internal  damping  and  elastic  forces  can  be  expressed  In 
terms  of  displacement  derivatives  and  material  constants. 


F'E(t)  -  Kx(t)  -  Cx(t)  =  Mx(t) 
where:  k  »  elastic  constant  c  =  damping 


By  definition  a  free  vibrating  body  Is  absent  of  external 
forces.  In  addition,  the  energy  loss  mechanism  for  dry 


I  1  .  Theor^ 


Finite  Elements 


k.  . 


Fg.(t)  -  FgCt)  -  Fc(t)  =  Mx(t) 
where: 


x(t)  = 
F^Ct)  * 


acceleration 
external  force 


Fc{t)  =  internal  elastic  forces 


Fc(t)  = 
M  = 


damping  forces 
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friction  and  hysteresis  losses  Is  not  well  understood  and 
as  Is  often  done  with  high  stiffness  low  mass  structures  the 
structural  damping  was  neglected  In  this  analysis.  Assuming 
displacements  are  harmonic,  equation  (2)  becomes  the 
fundamental  statement  of  the  vibration  problem. 


{K-u^M}0  =  0 
where: 

u  =  circular  frequency  (rad/ sec) 
D  =  displacement  vector 


(3) 


The  essence  of  the  FEM  Is  the  redefinition  of  a  continuous 
structure  In  terms  of  displacements  at  discrete  points  or 
nodes.  The  STAGS-Cl  code  Is  an  energy  based  code  using  the 
K1 rchoff-Love  hypothesis  [1].  For  a  conservative  system  In 
equilibrium  the  total  potential  energy  must  be  stationary 
and  the  first  variation  will  equal  zero.  The  total 
potential  energy  of  a  body  Is  the  strain  energy  of  the  body 
minus  the  work  done  on  the  body  by  external  forces.  Thus, 
the  total  potential  energy  for  a  free  vibrating  body  Is: 


n  =  u 


where 


^  J^olume  +  Vx  ' 


^yYxy)dVol 


(4) 
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The  kinematic  relations  In  STAGS  are  based  on  Sander's  shell 


equations  [1]  expressing  total  strain  as  a  linear  variation 
of  extenslonal  (membrane)  and  bending  (curvature)  strains  at 
the  midsurface.  The  midsurface  curvatures  being  composed  of 
components  of  rotation  about  coordinate  lines  and  the 
surface  normal . 


The  midsurface  strains  are  [2]: 

1/2  +  1/2  *2 

ey°  =  V  y  +  W/R  +  1/2  -  1/2  ^2 

where;  ^y,  and  ♦  are  the  components  of 

rotation  about  coordinate  lines  and  about 


the  normal  to  the  surface.  (Fig  6) 


(5) 

(6) 
(7) 


The  midsurface  curvatures  are: 

K  =  4 
X  ^x,x 


K  =  4 

y  ^y.y 

2K  =  2K 
xy  yx 


♦y.x  ♦x.y 


(8) 

(9) 

(10) 


and  In  terms  of  displacements: 


♦y  =  -W^y  +  V/R 

♦  * 


(11) 

(12) 


1/2  (U  „  -  V  ^) 


(13) 


Through  use  of  the  K 1 rchoff-Love  hypothesis  and  the 


orthoropic  constitutive  relations,  the  strain  In  any  layer 
of  the  laminate  can  be  expressed  In  terms  of  the  previously 
derived  midsurface  strains  and  curvatures 


where  [Q]j^  Is  the  transformed  stiffness  and  Is  composed  of: 


(15) 


In  which 


Qll 

=  QjjC^  +  2(Qj2  +  ^22®"' 

(  16) 

Qi2 

=  (.Q,l  +  Q22  -  4Q66)s2c2  Q,2(5^-"c^> 

(17) 

Q22 

=  QijS^  +  2(Qj2  +  2Qgg)s2c2  +  0^20^ 

(18) 

Qi6 

=  (Qii  -  Qi2  “2Q66)sc^  +  (Qj2  -  Q22 

+  2Qgg)s3c 

(19) 

7 


^26 


^66 


(Qll  -  Qi2 

(Ql 1  +  Q22 


-  2Q,^)s"‘c^  +  (Q 


<66 


12 


^22  2Qg6)sc' 


(20) 


2Qi2  -  2Q6g)s^c2  + 


<66 


(s^  +  c^) 


(21) 


where  s  =  sine  and  c  =  cos0  (Fig  1).  The  angle  8  being 
defined  as  the  angle  from  the  x  axis  of  the  body  coordinate 
system  to  the  1  axis  coincident  with  the  fibers. 


Q,  1  =  Ej/(  I-V21  '»i2^ 


^  “  ''12^2/^^  ”  ''12''21^  *  ''21^1^^*  ■  '<12''21^ 


^22  * 


12''21 


^66  =  ^12 


(22) 

(23) 

(24) 

(25) 


,  • .  •  . 


V.  - -1 


By  integrating  through  the  thickness  of  the  laminate  and 
remembering  midsurface  strains  and  curvatures  are  not  a 
function  of  z  [3]  the  resultant  forces  and  moments  (per  unit 
length)  can  be  expressed  as; 


(See  Figure  2  for  sign  notation  for  the  forces  and  moments 


and  Figure  3  for  the  geometry  of  the  laminate.) 

The  strain  energy  of  a  panel  can  then  be  expressed  In  terms 
of  moment  and  forces. 


Or  substituting 


equations  26-27: 


IJ 

«IJ 

1J 

°1J 

dArea  (32) 


The  finite  element  Is  constructed  from  assumed  displacements 
fields  that  serve  as  Interpolation  formulas.  The  exact  and 
assumed  displacement  (Interpolated  curve)  match  at  the  nodes 
but  may  vary  elsewhere.  This  however  does  not  Imply  nodal 
values  are  exact  [4],  The  plate  elements  In  STAGS  Cl  use  an 
Hermit  I  an  Interpolation  In  which  the  ordinates  and  slopes  at 
the  end  po’nts  of  adjacent  elements  are  forced  to  match. 

This  operation  results  In  an  evaluation  of  the  Interpolation 
functions  at  the  node  ordinates.  This  results  In  the 
formulation  of  shape  functions,  which  define 
displacements  within  an  element  when  the  Ith  element  degree 


of  freedom  has  untt  value  and  all  other  degrees  of  freedom 


are  zero.  The  strains  In  equation  (32)  can  be  expressed  In 
terms  of  nodal  displacements  using  the  strain-displacement 
matr lx  [G] . 


[G]  {d} 


(33) 


where  [G]  Isa  matrix  populated  with  first  order  derivatives 
for  Inplane  strains  and  second  order  derivatives  for  bending 
strains.  finally,  substituting  equation  (33): 


n  a  1/2 

f  £d}^ 

[G]"^ 

’*u 

«1J 

B.J 

°1J_ 

el ementa 1 

stiffness  matrix  Is  then 

[kij]  = 

f 

*1  j 

®1J 

[G] 

> 

°1J_ 

[G]{d}  dA 


(34) 


(35) 


and  Is  superimposed  with  the  other  elemental  stiffness 
matrices  to  form  the  global  stiffness  matrix  for  the  panel. 

The  mass  matrix  used  In  the  STAGS  Cl  code  1 s  of  a 
lumped  diagonal  form  [13].  Conceptually,  lumping  Implies 
that  translational  motion  of  a  node  Is  shared  evenly  by  the 
adjacent  nodes.  The  mass  can  be  thought  of  as  equivalent 


1  umps  at  the  nodes  and  the  1 umps  have  no  rotary  1 nert 1  a . 


The  implication  of  zero  rotary  inertia  is  that  zeros  are 
located  on  the  diagonal  corresponding  to  rotational  degrees 
of  freedom.  Lumped  mass  matrices  are  positive  semi -def i n i te 
if  they  have  zeros  in  the  diagonal.  The  implications  being 
that  convergence  from  above  (computed  natural  frequencies 
greater  than  exact)  cannot  be  assured.  In  fact,  natural 
frequencies  using  lumped  mass  matrices  usually  yield  natural 
frequencies  lower  than  exact  [4]. 

By  defining  the  structural  and  mass  properties  in  terms 
of  displacements,  the  solution  of  equation  (3)  becomes  that 
of  an  eigenvalue  problem.  For  nontrivial  solutions  the 
determ i nant 

[K,j]  -  \  [M,j]  =  0  (36) 

•p 

where  X  =  w 

The  determinant  results  in  the  characteristic  equations 
whose  roots  are  known  as  eigenvalues  (X).  Substitution  of 
the  eigenvalue  Into  equation  (36)  results  In  a  matrix  of 
coefficients  defining  the  eigenvector  or  natural  mode.  The 
eigenvector  describes  the  relative  displacements  of  the 
nodes  in  n-dimensional  space.  It  Is  Important  to  note  that 


the  eigenvector  only  defines  relative  displacements.  Hence, 


the  shape  of  the  natural  modes  is  unique,  but  not  the 
ampl 1 tude  [ 7] . 

Holographic  1 nterf erometv 

A  hologram  is  formed  when  two  sets  of  spatially 
coherent  light  waves  are  made  to  interfere  on  a  photographic 
medium.  The  light  wave  originating  at  a  single  light  source 
is  split  between  two  paths;  one  illuminating  the  object 
being  recorded,  the  other  serving  as  a  reference  to  the 
original  coherent  wavefront.  The  wave  front  reflected  from 
the  recorded  object  is  irregular  and  very  complicated.  On 
the  reference  path  the  spatial  coherence  is  retained.  The 
two  paths  intersect  again  at  the  photographic  plate  (Fig  4) 
creating  a  complex  interference  pattern.  The  interference 
is  of  two  types,  constructive  and  destructive.  When  the  two 
wavefronts  meet  at  a  point  in  phase,  the  total  amplitude 
increases  constructively  as  the  sum  of  the  Individual 
amplitudes.  Simlliarly,  when  the  two  wavefronts  are  not  in 
phase  a  destructive  or  diminishing  effect  occurs  in  the 
amplitude.  This  interference  pattern,  known  as  a  field,  is 
recorded  on  the  photographic  plate  during  the  exposure.  The 
points  where  constructive  Interference  occurs  form  parallel 


dark  fringes  on  the  order  of  5  x  10* 


[5]  and  are  not 


visible  to  the  eye.  After  the  plate  Is  developed  and  fixed 
these  fringes  become  the  hologram  and  act  much  like  a  wave 
grating  [16].  When  this  hologram  Is  Illuminated  by  a 
coherent  light  source  Identical  to  the  reference  beam,  the 
fringes  diffract  the  light  Into  wavefronts  essentially 
Indistinguishable  from  the  waves  which  were  reflected  from 
the  object.  (Fig.  5)  The  diffracted  waves  produce  all  the 
optical  characteristics  of  phase  and  amplitude  of  the 
original  waves  [5].  This  Image  differs  dramatically  from  a 
normal  photograph  which  Is  essentially  a  spatial  recording 
of  light  Intensity  on  the  film. 


Conduct  Inq  .-1  brat  Ion  analysis  with  holography  requires 
a  time  ave- aging  technique.  The  panel  Is  excited  (In  this 
case  acoustically)  at  a  natural  frequency  and  the  object 


assumes  a  modal  vibration  pattern.  The  photographic  plate 
records  all  fields  that  existed  during  the  exposure  time  In 
proportion  to  the  fraction  of  time  during  which  the 
wavefront  existed.  Even  though  the  fields  are  reconstructed 


mwaviv 


simultaneously  by  the  reference  beam  [5],  Thus,  when  viewed 
by  an  observer,  a  pattern  of  broad  fringes,  much  greater  In 
size  and  clearly  visible  to  the  eye,  appear  on  the  object 
Image.  Due  to  the  Interference  between  the  fields,  these 
broad  fringe  lines  connected  points  of  constant  deflection. 
It  should  be  noted  at  this  point  that  only  the  out  of  plane 
displacements  are  recorded  In  this  manner  and  care  must  be 
taken  to  avoid  Interpreting  slmlllar  deflection  patterns  on 
a  single  body  as  necessarily  Indicating  the  same  mode.  Tw 
modes  may  have  slmlllar  out-of-plane  displacements  yet  have 
In-plane  displacements  which  are  dramatically  different 


MODELLING 


111. 


Panel  Properties. 

The  graphite/epoxy  shell  study  In  this  thesis  was 
cylindrical  In  shape  with  a  height  of  16  Inches,  a  chord 
length  of  16  Inches,  and  a  radius  of  curvature  of  12  Inches. 
The  shell  dimensions  and  sign  convention  used  In  the  STAGS 
program  are  shown  In  (Fig  6).  The  panels  were  eight  quasl- 
Isotroplc  layups.  [0,  -45,  +45,  90]s,  of  Hercules  3501 /AS 
graphite  epoxy.  Each  ply  Is  approximately  .005  Inches  thick 
for  a  total  laminate  thickness  of  .04  Inches.  The  cure 
cycle  can  be  found  In  Appendix  B.  The  fabrication  was 
conducted  by  personnel  at  the  Air  force  Flight  Dynamics 
Laboratory  as  were  the  material  property  tests.  These  tests 
determined  the  material  properties  of  the  material  as; 


^1 

= 

18.84  X  10®  PSI 

^2 

= 

14.68  X  10®  PSI 

*^12 

= 

9.099  X  10®  PSI 

-2 

'*21 

S 

2.1813  X  10  ^ 

P 

s 

.055  lb/1n^ 

12 


=  .2799 


These  parameters  were  an  Input  to  the  STAGS-Cl  code.  It 
should  be  noted  that  the  User's  Manual  for  STAGS  calls  for 
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thickness  =  8  plies  @  0.005”  =  0 
radius  =  12” 

width  =  chord  length  =  16" 
length  =  16" 

structural  coordinate  directions 
ply  orientation 
displacements 
rotat I ons 

F I gure  6 :  Pane 1  Notat 1 on 


Vj2f  but  this  notation  is  contrary  to  that  found  in  most 
references  for  composite  materials.  The  property  expected 
by  STAGS  ^21  conventional  notation 

'*21  ~ 

Where  ’s  defined  as  the  strain  normal  to  the  fiber 
direction  and  Cj  is  defined  as  strain  in  the  fiber 
dl rect I  on . 

E 1 ement  Se 1 ect i on 

Curved  shell  elements  require  large  amounts  of  computer 
time  for  formulation  of  the  stiffness  matrix.  Also  the  use 
of  a  curved  element  can  Induce  strain  under  rigid  body 
motion  [0].  Thust  the  approach  used  by  the  STAGS  code  Is  to 
use  flat  plate  elements.  These  QUAF  elements  incorporate 
both  in-plane  and  bending  degrees  of  freedom.  Modelling  a 
shell  with  flat  plate  elements  can  result  In  element 
nonconformity  between  adjacent  elements.  The  flat  elements 
meet  at  angle  and  a  normal  rotation  degree  of  freedom  Is 
added  at  the  corner  nodes  to  enforce  rotational 
compatibility  of  element  boundaries  [Fig.  7]. 


Th ! s  Implies: 


(Byj  -  -  (B^i  ■•■  ^■22^s\T\{.<i/Z)  =  0 

-  ^^2'^cos{<x/Z)  -  (Byj  +  By2)5in(o/2)  =  0  (38) 

where : 

By  =  rotation  about  y  coordinate  line 

B^  =  rotation  about  z  coordinate  line 

a  =  angle  between  elements 

The  STAG-Cl  code  defines  a  small  limit  and  if  a  <  qq,  the 
rotation  B^  is  ignored  and  as  an  approximation  the 
conformity  constraint  becomes  By^  =  By2  [8]» 

A  more  serious  problem  is  caused  by  interelement 
displacement  compatibility.  The  assumed  displacement  fields 
for  the  in-plane  displacementst  u  and  v,  contain  first  order 
derivatives  while  the  transverse  displacement  field  for  w 
contains  second  order  derivatives.  Therefore,  w  is  usually 
represented  by  higher  order  ploynomials  than  u  and  v. 

Because  the  elements  meet  at  an  angle,  compatibility  between 
the  boundaries  of  adjacent  elements  require  that: 

(Vj  -  V2)cos(q/2)  -  (Wj  +  W2)sin(a/2)  ■  0 
(w,  -  W2)co5(a/2)  +  (v.  +  V2)sin(a/2)  *  0 


(39) 


of  these  displacement  derivatives  as  degrees  of  freedom 
result  In  an  In-plane  displacement  field  that  Is  cubic  for 
displacements  normal  to  element  boundaries.  The 
displacements  tangential  to  the  element  boundaries  are  made 
quadratic  by  raising  the  order  of  u  In  the  x  direction  and  v 
In  the  y  direction  through  the  addition  of  tangential 
degrees  of  freedom  at  midside  nodes.  The  midside  nodes 
allow  linear  strain  tangential  to  the  boundaries  and 
alleviate  the  slow  convergence  found  with  constant  strain 
elements  [8].  The  complete  set  of  shape  functions  In 
Isoparametric  form  can  be  found  In  [8].  The  STAGS-Cl  QUAF 
411  Incorporates  these  7  degrees  of  freedom  at  the  corner 
nodes  and  1  degree  of  freedom  at  each  of  the  midside  nodes 
for  a  total  of  32.  (Fig  8)  A  simpler  QUAF  410  element  (Fig 
9)  Is  available  that  eliminates  the  side  nodes  and  uses  only 
an  average  normal  rotation  at  corner  nodes.  This  reduces 
the  element  to  constant  strain  capability  tangential  to  the 
boundaries  and  suppresses  shear  strain  at  the  corner  nodes. 
Consequently  one  would  expect  faster  convergence  with  fewer 
elements  when  modelling  curved  surfaces  with  the  QUAF  411 
e 1 ement . 
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IV.  FINITE  ELEMENT  ANALYSIS 


Defining  the  Element  Grid 

At  the  beginning  of  this  effort.  It  was  felt  that  the 
QUAF  411  plate  element  would  be  required  to  model  the  curved 
surface  of  the  panels.  This  decision  was  based  on  the 
normal  degrees  of  freedom  added  to  this  element  to  Insure 
compatibility  between  element  boundaries  subject  to  shear. 
The  Inter element  compatibility  requirement  and  higher  order 
normal  strain  (linear  vs  constant)  function  were  predicted 
to  provide  faster  convergence  with  fewer  elements.  A 
potential  liability  of  this  choice  was  that  the  32  degrees 
of  freedom  and  nine  point  numerical  Integration  utilized  by 
the  411  element  versus  the  24  degrees  of  freedom  and  5  point 
Integration  of  the  410  [8]  element  would  require  greater 
computer  resources  for  the  same  grid  size.  This  was  of 
particular  concern  due  to  the  desire  to  utilize  a  digital 
VAX  11/780  computer  operating  under  the  VMS  operating  system 
for  analysis. 

A  convergence  study  was  conducted  using  a  32  x  13,  19  x 
19,  and  25  x  25  square  element  grid  for  the  solid  shell. 

This  results  In  a  constant  1.0,  0.75,  and  0.5,  Inch  grid 
line  spacing  respectively.  The  results  are  tabulated  below 


for  the  fundamental  frequency. 


4 1 1  Computat  f  on  T f  me  410  Computat Ion  Time 


GRID 

FREO 

CPU  sec 

FREO 

CPU  sec 

13  X 

13 

364  Hz 

698 

365  Hz 

439 

19  X 

19 

460  Hz 

2355 

461  Hz 

1413 

25  X 

25 

504  Hz 

5844 

506  Hz 

3019 

E 1 ement  Compar 1  son 
Table  1 

It  Is  apparent  that  the  anticipated  benefits  of  the  411 
element  did  not  materialize.  The  computed  fundamental 
frequencies  varied  by  less  than  0.5  percent.  The  other 
result  of  Interest  lies  In  the  difference  In  computation 
time.  The  computation  time  for  the  411  element  Is  37 
percent  greater  for  the  13  x  13  grid  and  48  percent  greater 
for  the  25  x  25  grid.  The  410  element  appears  to  provide 
virtually  Identical  results  with  a  substantial  cost  savings. 
The  failure  of  the  411  element  to  exhibit  faster  convergence 
suggests  that  for  the  c 1 amped-c 1 amped  boundary  conditions 
the  higher  order  displacement  fields  are  Ineffective.  This 
can  be  explained  by  reviewing  the  displacement  fields  of  the 
additional  degrees  of  freedom  found  In  the  411  element.  The 
corner  node  rotation  Is  added  to  the  element  as  derivatives 


of  Inplane  displacements  with  respect  to  the  normal  ordinate 


i.e.,  U,  and  V,  .  The  midside  nodes  also  add  only  to  the 
y  ^ 

order  of  the  Inplane  displacement  fields.  Considering  the 
boundary  conditions,  the  small  vibration  problem  is 
primarily  one  of  out-of-plane  displacements.  Thus,  the 
additions  to  the  order  of  the  inplane  displacements  add 
little  to  the  solution  of  this  particular  problem.  The  QUAF 
410  element  was  Judged  best  for  this  problem. 

To  insure  the  grid  size  had  been  reduced  to  a  point 
where  the  tradeoff  between  accuracy  and  computation  time  had 
reached  an  acceptable  level,  a  computer  run  utilizing  a  31  x 
31  square  grid  was  submitted.  This  run  failed  to  converge 
to  acceptable  error  limits  for  the  fifth  eigenvalue  before 
the  al loted  disk  space  on  the  VAX  11/780  was  exceeded. 
However,  the  fundamental  eigenvalue  had  converged  and  was 
determined  to  be  529  Hz.  This  is  a  change  of  approximately 
5  percent,  and  the  25  x  25  grid  was  established  as  the  best 
balance  between  accuracy  and  computer  resources. 

As  indicated  in  the  previous  discussion,  a  decreasing 
grid  mesh  results  in  greater  computed  natural  frequencies. 
Since  the  finite  element  model  is  an  approximation  utilizing 
a  finite  number  of  degrees  of  freedom  to  model  an 
essentially  infinite  number  found  In  the  actual  structure, 
the  model  should  always  be  expected  to  possess  greater 


stiffness.  Given  the  same  mass,  the  model  would  then  be 


expected  to  yield  higher  natural  frequencies  with  Increasing 
grid  spacing  (fewer  degrees  of  freedom).  Conversely,  as  the 
spacing  Is  refined  the  computed  frequencies  should  decrease 
or  converge  from  above  to  the  exact  value.  The  behavior 
Indicated  In  this  analysis  would  then  Imply  some  error  In 
the  model.  As  mentioned  In  the  theory  discussion,  the 
lumped  mass  matrix  frequently  causes  the  results  to  converge 
from  the  bottom.  This  behavior  was  also  observed  In  an 
evaluation  of  the  STAGS  Cl  code  by  the  Office  of  Naval 
Research  [15].  This  study  noted  that  while  the  convergence 
problem  Indicated  a  failure  to  satisfy  some  equations  of 
mechanics,  the  results  were  not  adversely  affected. 

Panel  Analysis 

Having  determined  the  grid  size  to  be  used  for  each 
panel,  an  Input  file  for  each  of  the  panels  was  constructed. 
The  cutouts  were  entered  using  the  IRREG  parameter  on  the  N- 
1  Input  card.  Then  the  N-8  Input  card  allows  the  user  to 
define  rectangular  cutouts  In  terms  of  grid  coordinates  [6]. 
The  2x2  Inch  square  being  defined  row  11-15,  column  11-15; 
the  2x4  cutout  row  11-15,  column  9-17;  the  4x4  cutout 
row  9-17,  column  9-17.  To  facilitate  the  plotting  of  the 


K  .i 


eigenvectors,  the  z  axis  of  the  global  reference  system  was 
chosen  to  intercept  the  midpoint  of  the  chord  of  the  panel. 
(Fig  10). 


This  placement  allowed  the  displacements  expressed  In  terms 
of  the  local  curvilinear  system  to  be  easily  rewritten  based 
on  the  global  system.  This  conversion  was  conducted  using 
the  expression 

w  cos{y)  +  V  s1n(y)  =  W  (40) 

Where  y  =  local  coordinate  (degrees) 
w  =  local  normal  displacement 
V  =  local  transverse  displacement 
W  =  global  displacement 

The  transformed  displacements  could  then  be  plotted  using  a 
contour  plot  algorithm  and  provide  the  same  perspective  of 
mode  shape  found  In  the  holograms.  The  global  system  was 
fixed  by  Inputing  the  appropriate  global  value  for  the 
extreme  corner  points  as  required  by  STAGS.  Remembering  the 
Y  coordinate  line  Is  expressed  In  terms  of  degrees  the 
following  comparison  Is  Included. 


Corner  Points 
C 


Local  System 


Global  System 


1 

1 

0.0 

-28.65 

0.0 

-5.75 

10.53 

1 

25 

0.0 

28.65 

0.0 

5.75 

10.53 

25 

1 

12.0 

-28.65 

12.0 

-5.75 

10.53 

25 

25 

12.0 

28.65 

12.0 

5.75 

10.53 

Global  System  Definition 
Table  2 

A  sample  Input  file  for  a  panel  with  and  without  a  cutout  Is 
1 nc 1 uded  f n  Append lx  A . 

Numerical  Solution 

The  STAGS-Cl  code  utilizes  an  Inverse  F>ower  method  to 
evaluate  the  roots  of  equation  (36).  This  method  [4] 
rewrites  equation  (36)  as: 

[K]  -  Xq[M]  (D)  =  [M]  {D')  j  (41) 


tD')p  =  Dp/Cp 


where 


X  =  shift  point 
o 

Cp  =  largest  magnitude  of  vector  {D}p. 


The  algorithm  assumes  a  vector  solution  (0'}^  (p=l)  to 

solve  equation  (41)  for  {D}j.  {0} j  is  then  scaled  using 

equation  (42)  to  get  {D')j  which  in  turn  is  reiterated  in 
equation  (41).  This  iterative  method  converges  to  the 
eigenvalue  closest  to  the  shift  point  Aq.  Since  in  this 
analysis  the  lowest  five  eigenvalues  were  desired,  the  shift 
point  (SHIFT  on  input  card  D-3)  was  chosen  as  zero.  The 
iteration  values  can  be  included  as  part  of  the  program 
output  (1  PR I  NT  on  D-2  card).  The  eigenvalues  are  then  used 
to  evaluate  equation  (36)  to  determine  the  eigenvectors. 

The  program  lists  the  displacements  in  tabular  form  at  each 
grid  point.  These  displacements  were  edited  in  a  short 
Fortran  routine  using  equation  (42)  and  the  displacement 
values  rewritten  in  the  data  field  format  required  by  the 
contour  plot  routine.  The  plot  routine  was  then  exercised 
to  provide  a  visual  outline  of  the  mode  shape  at  each 


natural  frequency. 


S  S  V 


V.  HOLOGRAPHIC  ANALYSIS 


i! 


I 

■ 


The  Test  Fixture 

In  order  to  conduct  the  holographic  analysis,  a  test 
fixture  capable  of  supplying  clamped  boundary  conditions  on 
both  the  flat  and  curved  edges  of  the  cylindrical  shell  was 
needed.  In  addition,  the  test  fixture  had  to  remain  stable 
over  the  frequency  range  of  Interest  and  under  temperature 
changes  that  might  occur  during  testing.  The  extreme 
sensitivety  of  the  holographic  method  Is  quite  capable  of 
recording  rigid  body  motion  of  the  test  fixture  due  to 
thermal  effects  or  vibration.  This  rigid  body  motion  would 
be  superimposed  on  the  panel  motion  and  would  be  Impossible 
to  distinguish.  Steel  was  chosen  over  aium1n*im  due  to  Its 
lower  coefficient  of  thermal  expansion.  gi  eater  stiffness, 
and  mass.  The  fixture  was  designed  and  constructed  to 
provide  a  two  Inch  clamping  surface  on  eac.  ige .  The 
vertical  clamps  were  constructed  (f’g  li)  *  \2  inch 

sections  of  4  x  4  Inch  bar  stock.  A  £  I /0  *  I  1/2  path  was 
milled  from  one  face  of  each  bar  to  form  a  u  shaped 
structure.  Two  rows  of  quarter  Inch  holes  were  drilled  and 
tapped  along  each  side  for  clamping  bolts.  These  bolts 
tightened  against  1/2  x  2  x  12  inch  flats  placed  on  each 
side  of  the  graph  I te/epoxy  panel.  When  the  bolts  on  each 
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stde  of  the  clamp  were  tightened  the  panel  was  securely 
clamped  between  the  flats.  To  prevent  lateral  motion 
smaller  U  shaped  clamps  were  Installed  on  the  open  face  of 
the  vertical  clamps.  A  curved  1  1/4  Inch  wide  by  2  1/8  deep 
path  with  an  outer  radius  of  twelve  Inches  was  milled  from  a 
22  Inch  wide  by  10  Inch  deep  by  3  Inch  thick  steel  stack. 

The  twelve  Inch  outer  radius  matched  that  of  the  panels.  A 
set  of  9  curved  blocks,  1  Inch  thick,  2  Inches  In  height, 
and  1  3/4  Inches  In  chord  were  milled  to  match  the  Inner 
radius  of  the  panels.  Great  difficulty  In  accurately 
machining  these  blocks  lead  to  a  decision  to  add  a  thin 
(40  mils)  ductile  aluminum  bushing  between  the  outer  radius 
of  the  block  sets  and  the  panel.  It  was  felt  that  this 
bushing  would  evenly  distribute  clamping  forces  between 
clamping  blocks.  Transverse  motion  here  was  prevented  by 
milling  two  lateral  grooves  to  match  the  contact  points  of 
the  clamping  bolts.  The  paths  for  these  bolts  were  drilled 
and  tapped  radially  on  the  Inner  radius  of  the  horizontal 
flats.  Four  1/2  Inch  rods  were  added  to  provide  structural 
rigidity  and  act  as  assembly  aids  when  changing  panels. 

When  assembling  the  fixture,  the  two  vertical  clamps  and 
bottom  section  were  Joined  and  the  panel  Inserted  between 
clamping  surfaces.  The  top  flat  was  then  Installed  and  the 
radial  clamping  blocks  were  Inserted  from  the  concave  side. 


Each  clamping  bolt  was  tightened  to  6  ft- lbs.  This  fixture 
appeared  to  provide  a  very  good  approximation  of  the  clamped 
boundary  conditions  that  were  desired.  In  none  of  the 
holography  experiments  were  fringe  lines  ever  observed  to 
extend  Into  clamp.  Such  a  fringe  line  would  have  Indicated 
a  line  of  displacement  extending  Into  the  clamp  and  thus 
poor  c 1  amp  1 ng . 


Experimental  Procedure 

The  equipment  used  In  the  experiment  Is  listed  below: 


Function  Generator: 
Frequency  Counter: 

Audio  Ampl 1 f 1 er : 

Fotonic  Displacement  Meter: 

Osci 1 loscope: 

50mV  Hellum/Neon  Laser: 

Acoustic  Drivers: 


H  P  Model  331  OB 

H  P  Model  5316A 

Bogen  Model  HTA125 

Mechanical 

Technology  KD  320 

Tektronix  7603 

Spectre  Physics 

(Model  125) 

Atlas  Sound  60  Watt 


A  block  diagram  of  the  equipment  configuration  Is  shown  In 
Fig  (12).  The  opt  1 ca 1  set  up  I s  shown  In  Fig  (13).  In  order 
Illuminate  the  entire  area  of  the  panel  and  photographic 


plate,  20X  microscope  objectives  were  employed  to  expand  the 
beam.  Pinhole  spatial  filters  are  used  In  series  with  the 
objectives  to  eliminate  refractive  fringes  due  to  dust 
particles  on  the  mirrors. 

Real  time  holographic  Interferometry  was  used  to 
determine  the  natural  frequencies  and  mode  shapes  of  the 
four  panels.  The  procedure  Involves  making  a  hologram  of 
the  stationary  panel.  The  hologram  Is  then  placed  In  the 
plate  holder  In  the  position  In  which  It  was  formed.  When 
this  hologram  and  the  object  are  Illuminated,  their  Images 
are  coincident.  When  the  panel  Is  excited  at  resonant 
frequency  by  the  acoustic  drivers.  It  assumes  a 
corresponding  mode  shape.  The  two  Images  are  no  longer 
coincident.  Thus,  broad  fringe  bands  are  formed  on  the 
Image  Into  a  pattern  which  resembles  that  of  the  associated 
mode  shape.  The  procedure  then  becomes  one  of  sweeping 
through  the  frequency  range  until  a  pattern  appears  on  the 
Image.  To  preclude  the  possibility  that  the  pattern  Induced 
by  the  Input  signal  of  the  acoustic  drivers  Is  not  actually 
a  response  to  a  multiple  (harmonic)  of  that  frequency,  an 
optical  displacement  meter  was  employed  to  determine  panel 
response.  By  comparing  the  output  signal  of  the  optical 
displacement  sensor  to  that  of  the  acoustic  signal,  the 
Ltssajous  pattern  that  results  on  the  oscilloscope  was  used 


to  confirm  that  the  panel  was  actually  vibrating  at  the 
acoustic  frequency.  For  example,  a  harmonic  response  to  a 
60  Hz  Input  signal  may  actually  be  occurring  at  480  Hz.  To 
the  viewer  of  the  hologram,  the  displacement  patterns  are 
Identical  but  the  Lissajous  pattern  would  actually  display  8 
lobes.  The  displacement  sensor  had  two  advantages  over  an 
accelerometer  often  employed  In  this  instance.  The  optical 
probe  was  easily  moved  about  the  face  of  the  panel  without 
having  to  break  any  bonding  glue.  The  second  advantage, 
particularly  for  vibration  analaysis.  Is  that  the  optical 
technique  does  not  add  mass  to  the  structure.  The  added 
mass  associated  with  accelerometers  actually  charges  the 
response  of  the  structure  during  the  analysis. 

In  this  experiment,  two  Identical  acoustic  drivers  with 
exponentially  convergent  cones  were  used  to  excite  the 
panels.  Two  horns,  driven  180  degrees  out  of  phase,  were 
requ I  red  to  exc 1 te  mode  shapes  that  were  ant  1 symmetr 1 c  about 
the  vertical  line  of  symmetry  of  the  curved  panel.  This 
push-pul  1  excitation  was  particularly  necessary  when  a 
symmetric  and  antisymmetric  mode  shapes  occurred  at 
relatively  close  frequencies.  It  should  be  noted  that  the 
horns  used  In  the  experiment  were  approximately  6  Inches  In 
diameter  and  thus,  could  not  be  placed  on  the  same 
horizontal  plane  and  still  excite  the  concave  surface  of  the 


panels.  The  horns  had  to  be  placed  at  opposing  corners  and 
thus  were  not  capable  of  a  true  antisymmetric  excitation. 
This  proved  to  be  of  minor  consequence  as  discussed  in  the 
next  section. 

After  using  the  real  time  technique  and  optical 
displacement  sensor  to  determine  natural  frequencies,  a  time 
average  hologram  was  constructed  to  permanently  record  each 
mode  shape.  The  mode  shapes  and  frequencies  were  then 
compared  with  numerical  results. 


VI.  DISCUSSION  AND  RESULTS 


This  chapter  Is  divided  Into  four  sections.  It  Is 
Important  to  note  that  the  comparrison  offered  here  Is 
between  experiment  and  numerical  approximation.  Therefore, 
errors  both  positive  and  negative  can  be  expected.  Had  the 
comparison  been  a  closed  form  or  exact  solution,  the  grid 
spacing  would  be  expected  to  predict  values  less  than  exact 
under  al 1  conditions.  The  experimental  test  fixture  can 
only  approximate  clamped  conditions  that  are  modelled 
exactly  In  the  finite  element  solution.  Thus,  both 
solutions  studied  here  are  approximations  of  the  exact 
solution  and  no  trend  In  error  should  be  expected.  The 
results  for  each  of  the  four  panels  are  discussed  In  order 
of  Increasing  cutout  size. 

The  Solid  Panel 

The  solid  panel.  Its  shape  being  the  simplest,  was 
chosen  first  for  both  the  numerical  and  experimental 
analysis.  Since  the  lowest  five  natural  frequencies  were 
sought,  the  frequency  range  of  Investigation  started  at  zero 
and  progressed  upward  until  all  five  frequencies  were  found. 
Excitation  was  provided  by  a  single  horn  aimed  at  the  center 


of  the  panel.  The  soHd  panel  was  particularly  susceptible 
to  displaying  harmonics.  Using  the  real  time  holographic 
technique,  an  Identical  modal  displacement  pattern  was  seen 
to  repeat  Itself  between  0  and  534  Hz.  Through  use  of  the 
optical  displacement  meter.  It  was  determined  that  each  of 
the  Identical  displacement  patterns  below  534  Hz  were 
harmonic  responses.  This  technique  established  that  the 
true  natural  frequency  was  Indeed  534  Hz.  It  was  possible 
to  visibly  observe  this  process  by  scanning  the  frequency 
range  at  reduced  sound  levels  from  the  acoustic  horn.  The 
harmonic  responses  would  not  appear  at  sound  levels  that 
would  easily  excite  the  true  mode.  The  existence  of 
harmonics  Is  the  result  of  exciting  the  panel  with  point 
sources  such  as  the  acoustic  horns  equipped  with  the 
exponential  horns. 

Scanning  the  frequency  upward.  It  became  evident  that 
only  those  natural  frequencies  with  mode  shapes  symmetric 
about  the  vertical  midplane  were  being  excited.  A  second 
horn,  with  Input  leads  switched  to  produce  a  response  180 
degrees  from  the  other,  was  added  to  produce  an 
antisymmetric  excitation.  Scanning  again  from  zero,  a  modal 
response  antisymmetric  about  the  vertical  midplane  was  found 
at  525  Hz.  It  was  also  found  that  the  symmetric  mode  at  534 
Hz  could  also  be  excited  with  the  dual  horn  arrangement. 


The  small  difference  In  these  two  natural  frequencies,  1.7 
percent,  resulted  In  smearing  of  mode  shapes  In  the  hologram 
as  the  frequency  gap  was  traversed.  This  can  easily  be 
explained  due  to  the  fact  that,  while  the  natural 
frequencies  are  discrete  values  mathemat 1 ca 11 y ,  the  panel 
would  respond  simlllarly  to  a  small  range  of  frequencies 
about  that  point.  For  these  two  frequencies,  the  ranges 
were  observed  to  overlap.  The  frequency  scanning  was 
continued  In  this  manner  until  the  first  five  modes  were 
found.  A  general  observation  can  be  made  at  this  point  that 
those  mode  shapes  exhibiting  symmetry  about  the  vertical 
midplane  were  always  easier  to  excite. 

SOLID  PANEL 

Mode  Frequency  (Hz)  Percent 

Number  STAGS  Experiment  Error 


The  mode  shapes  determined  by  the  finite  element  method  and 
holographic  Interferometry  are  presented  In  Figures  14-18. 
The  fringes  represent  lines  of  constant  displacement  and  can 
be  Interpreted  much  like  a  topographical  map.  The  white 
areas  which  represent  nodal  lines  or  zones  of  zero 
displacement  can  be  deduced  from  the  fact  that  displacement 
1 s  zero  at  the  edges . 

As  a  comparison  of  the  photographs  and  contour  plots 
shows,  the  mode  shape  predictions  for  the  solid  panel  by  the 
STAGS-Cl  program  are  quite  accurate.  The  slight  tilting  of 
the  lines  of  constant  displacement  toward  the  upper  left  Is 
visible  In  both  the  holograms  and  contour  plots.  This 
effect  appears  to  be  a  bias  toward  the  +45  degree  direction 
which.  In  this  laminate,  corresponds  to  the  fiber  layers 
closest  to  the  midsurface  90  degree  lamina.  While  no 
explanation  for  this  behavior  Is  evident.  Its  occurrence  In 
both  cases  Indicates  It  Is  real  and  not  an  effect  of 
boundary  conditions  of  test  fixture.  Siml Marly,  the  error 
In  predicting  the  natural  frequency  Is  quite  low. 

The  2x2  Inch  Cutout 


As  with  all  the  panels,  the  procedure  employed  to  find 
natural  frequencies  was  Identical  to  that  used  for  the  solid 


Figure  16: 


panel.  The  first  and  second  natural  frequency  on  this  panel 
were  separated  by  approximately  7  percent  of  the  highest 
frequency.  The  smearing  exhibited  between  these  two  modes 
on  the  solid  panel  was  not  present  on  this  panel. 
Additionally,  the  harmonics  of  the  first  symmetric  mode 
(Mode  2)  were  as  evident  as  on  the  solid  panel.  The 
symmetric  mode  clearly  was  slmlllar  on  both  panels  with  the 
level  of  displacement  at  a  fixed  sound  level  significantly 
greater  than  the  antisymmetric  modes.  The  experimental  and 
numerical  results  are  tabulated  below. 

2x2  CUTOUT 


Mode 

Frequency  (Hz) 

Percent 

Number 

STAGS 

Experiment 

Error 

1 

505 

510 

-1.0 

2 

527 

546 

1 

UJ 

• 

Ul 

3 

669 

«« 

731 

-8.5 

4 

691 

705 

-2.0 

5 

712 

756 

-5.8 

Table  4 

•  • 


Order  Switch 


As  annotated,  the  order  (as  determined  by  the  mode  shape)  of 
modes  3  and  4  were  reversed  between  the  two  methods. 

However,  the  error  values  In  Table  2  are  computed  using  the 
frequencies  corresponding  to  Identical  eigenvectors  (Fig  21, 
22).  The  two  frequencies  of  mode  3  and  mode  4  differ  by 
less  than  4  percent  and  the  shift  may  be  attributed  to  less 
than  perfect  boundary  conditions  In  the  test  fixture. 
Problems  with  the  clamping  In  mode  3  and  4  are  visible  In 
the  holograms.  Note  that  the  regions  of  displacement  to  the 
right  of  the  cutout  do  not  descend  to  the  bottom.  Clamping 
appears  to  be  weak  on  the  bottom  clamp.  Despite  the  shift, 
the  computed  natural  frequencies  still  exhibit  good  accuracy 
for  all  five  modes.  The  eigenvectors  are  conpared  In 
Figures  19-23.  Again,  STAGS-Cl  proved  outstanding  In 
predicting  mode  shapes.  When  compared  with  the  mode  shapes 
of  the  solid  panel,  the  smallest  cutout  appears  to  have  very 
little  effect  for  the  first  four  modes.  Slmlllarly,  the 
shift  toward  the  upper  left  is  evident  in  some  modes. 
Although  mode  5  differs  dramatically  between  the  solid  and 
the  2x2  Inch  cutout  panels  the  computer  code  was  able  to 
accurately  predict  the  change. 


The  2x41 nch  Cutout 


The  STAGS-Cl  and  holographic  results  are  summarized  In 
Table  5. 


2x4  CUTOUT 


Mode 

Frequency  (Hz) 

Percent 

Number 

STAGS 

Exper 1 ment 

Error 

1 

454 

419 

+7.7 

2 

536 

463 

+  13.6 

3 

562 

579 

-2.9 

4 

631 

670 

-5.8 

5 

703 

714 

-1 .5 

Table  5 

This  panel  was  the  only  one  that  exhibited  significant  error 
for  the  first  two  modes.  Comparing  the  holograms  and 
contour  plots  does  not  Identify  a  region  where  boundary 
conditions  are  weak.  It  should  be  pointed  out  that  for  the 
first  two  panels  STAGS  Just  slightly  under-estimated  the 
natural  frequencies  of  the  first  modes  and  had  a  tendency  to 


increasingly  underestimate  them  at  the  higher  ones.  For  the 
2x4  cutout,  the  program  overestimated  the  lower  mode 
natural  frequencies.  The  trend  toward  an  increasing 
negative  percentage  error  produced  in  this  case  the  lowest 
errors  at  the  higher  modes.  Thus  it  appears  that  either 
STAGS  overestimated  the  panel  stiffness  or  the  clamped 
boundary  conditions  were  "soft"  on  al 1  edges. 

The  mode  shape  estimations  of  STAGS  were  particularly 
good  for  this  panel.  (Figs.  24-28)  The  tilting  of  the 
ant  inode  regions  noted  in  previous  panels  is  again  evident. 
Additionally,  the  bulging  of  the  ant  inode  region  about  the 
cutout  is  apparent  in  the  holograms  and  the  numerical 
results.  The  loss  of  stiffness  at  the  center  of  the  panel 
is  much  more  visible  for  this  panel  than  the  2x2  inch 
cutout.  The  deflected  regions  about  the  first  symmetric 
mode  (2)  extend  into  the  cutout.  The  symmetric  third  mode 
is  the  only  one  seen  in  any  of  the  penels,  while  mode  4  and 
5  appear  remarkably  similiar  in  shape  to  modes  3  and  5  of 
the  solid  and  2x2  inch  cutout  panels.  The  failure  of  this 
panel  to  follow  the  mode  shape  patterns  exhibited  by  the 
previous  two  was  an  indicator  of  the  behavior  found  in  the  4 
X  4  inch  cutout  panel. 


The  natural  frequencies  are  tabulated  In  Table  6.  It 
Is  apparent  that  the  STAGS  results  overestimate  the  first 
two  frequencies  (though  to  a  smaller  degree)  as  It  had  for 
the  2x4  Inch  cutout.  Again,  the  trend  toward  Increasing 
underestimation  of  higher  modes  Is  Indicated.  Despite  these 
observations,  the  overall  level  of  accuracy  Is  still  quite 
good.  It  should  be  noted  that  the  4x4  Inch  cutout 
represents  approximately  11  percent  of  the  surface  area. 

4x4  CUTOUT 

Mode  Frequency  (Hz)  Percent 

Number  STAGS  Experiment  Error 


Th?s  reduction  in  area  precipitated  a  9  percent  decrease  in 
natural  frequency  for  the  experimental  results.  However, 
the  first  mode  shape  for  this  cutout  is  symmetric  and  the 
second  mode  shape  is  antisymmetric.  This  is  a  reversal  of 
the  behavior  exhibited  in  the  first  three  panels.  The 
finite  element  model  was  again  accurate  in  predicting  this 
behavior. 

Locating  the  second  natural  frequency  for  this  cutout 
proved  quite  difficult.  When  scanning  through  the  frequency 
range  the  mode  shapes  1,  3,  4,  and  5,  as  predicted  by  STAGS, 
were  quickly  located  using  the  two  horn  excitation 
arrangement.  The  accuracy  of  STAGS  in  predicting  natural 
frequencies  and  in  particular  the  mode  shapes  of  the  first 
three  panels  indicated  that  perhaps  the  previous  procedures 
has  missed  a  mode.  In  order  to  increase  the  intensity  of 
the  excitation,  one  of  the  two  horns  was  moved  to  the  front 
surface.  After  rewiring  to  return  the  horns  to  in-phase 
motion,  the  experiment  was  again  repeated.  The  mode  shape 
seen  in  Fig  30  (note:  horn  in  photo)  is  the  response  seen  at 
510  Hz.  This  mode  shape  does  not  exhibit  the  correlation 
seen  at  modes  1,  3,  4,  and  5  but  was  definitely  present.  To 
preclude  this  response  as  being  a  harmonic,  the  optical 
probe  was  employed.  The  large  displacement  patterns  seen  in 
Fig  30  are  indeed  oscillating  at  510  Hz.  It  proved 


Impossible  to  eliminate  all  motion  along  the  free  edge  at 
this  frequency.  While  the  Information  achieved  through  use 
of  the  optical  probe  and  the  holograms  were  not  conclusive, 
the  transition  to  the  mode  shape  seen  In  the  photograph 
Indicated  an  attempt  by  the  panel  to  assume  a  mode  shape. 

It  Is  strongly  felt  that  a  natural  frequency  does  indeed 
exist  at  510  Hz  and  that  the  correct  combination  of  boundary 
conditions  and  excitation  was  never  found  to  Illicit  the 
response  predicted  by  STAGS.  All  other  mode  shapes  In  this 
panel  were  predicted  with  accuracy  found  In  other  panels. 
Mode  3  and  5  (Fig  31  and  33)  indicate  some  looseness  In  the 
bottom  clamp.  Overall,  the  results  for  this  panel  as  with 
the  others  shows  a  close  correlation  between  the  finite 
element  and  experimental  techniques. 
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Cutout  Effects 


As  seen  In  the  previous  discussion,  the  boundary 
conditions  generated  by  the  test  fixture  had  a  significant 
effect  on  the  mode  shapes.  As  seen  In  Figure  21  the  loss  of 
clamping  at  the  bottom  of  test  fixture  resulted  In  failure 
of  the  displacement  pattern  to  extend  to  the  bottom  of  the 
panel.  Given  the  narrow  range  of  frequencies  Involved  and 
the  closeness  of  some  of  the  natural  frequencies,  a  small 
deviation  from  the  exact  value  due  to  Inexact  boundary 
conditions  would  represent  a  significant  deviation.  This 
deviation  would  effectively  mask  any  trends  or  conclusions 
based  on  experimental  values.  Thus,  the  close  correlation 
of  the  experimentally  determined  data  (frequencies  and  mode 
shapes)  with  that  of  numerically  derived,  supported  the  use 
of  the  numerical  values  to  analyze  the  effects  of  cutouts. 
This  conclusion  Is  further  supported  In  that  an  analysis 
based  on  the  numerical  results  would  be  the  results  based  on 
mathematically  perfect  clamped  boundary  conditions.  Figure 
34  Is  a  plot  of  natural  frequencies  generated  by  STAGS  C-1 
for  each  of  the  four  panels.  Although  there  are  a  limited 
number  of  data  points  this  plot  clearly  shows  a  trend  toward 
lower  frequencies  (particularly  for  the  highest  three  modes) 
when  cutouts  are  added.  Figure  35  graphically  displays  the 


effect  of  mass  removal  on  each  mode.  Note  that  this  Isa 
comparison  of  the  modes  as  they  occurred  in  the  numerical 
sequence.  Mode  3  for  the  2x4  inch  cutout  is  not  the  same 
shape  as  Mode  3  of  the  2x2  inch  panel;  both  simply  are  the 
third  from  the  lowest  frequency  for  each  panel  respectively. 
If  in  a  simplified  analysis  the  panels  are  modelled  as  a 
single  degree  of  freedom  system,  then: 

“  =VT~ 


where  K  =  stiffness 

M  =  Mass. 

Then  the  apparent  decrease  in  frequency  (Fig  35)  for  the 
2x2  inch  cutouts  suggests  a  decrease  in  stiffness  versus 
mass.  Similiarly,  the  Increase  in  frequency  for  the  modes 
1,  3,  and  4  for  the  4x4  cutout  implies  mass  is  decreasing 
faster  than  stiffness.  Although  the  data  points  are  limited 
in  number,  the  trend  may  indicate  that  material  removal 
along  the  vertical  plane  of  symmetry  is  dominant  for  panel 
st i ffness . 

Referencing  Figures  47,48,54,55,61 ,62, 69, and  70  an 
apparent  shift  in  symmetry  occured  for  the  first  and  second 
mode  shapes  of  the  4x4  inch  cutout  panel.  While  this 
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result  was  unexpected,  it  was  also  observed  experimentally 


To  further  explore  this  event,  additional  finite  element 
runs  with  the  same  material  properties  were  completed  for  a 
3x3  and  3.6  x  3.6  Inch  cutout.  The  results  were  are  shown 
on  Fig.  35.  The  3x3  Inch  hole  retained  the  symmetry  seen 
In  the  smaller  cutouts  while  the  larger  3,6  x  3.6  Inch 
cutout  displayed  the  symmetry  shift  seen  In  the  4x4  Inch 
cutout.  Additional  runs  were  then  accomplished  In  order  to 
explore  the  contribution  of  material  properties  to  this 
apparer.t  shift  In  symmetry.  While  the  original  geometry  of 
the  panel  was  retained,  the  orthotropic  material  was 
replaced  with  an  Isotropic  material  with  a  modulus  of 
elasticity  equivalent  to  Ej.  The  symmetry  shift  was  again 
evident  for  the  4x4  Inch  cutout.  Thus,  If  the  symmetry 
shift  Is  Indeed  real  It  appears  to  be  a  function  of  the 
panel  radius  and  cutout  size. 

In  conclusion,  the  effects  of  the  Interior  cutouts 
tends  to  create  a  moderate  decrease  In  natural  frequency  for 
this  panel  geometry  and  material  properties.  The  geometry 
of  the  problem  analyzed  In  this  case  appears  to 
significantly  alter  the  mode  shapes  of  the  panels.  This 
observation  Is  contrary  to  that  found  for  flat  plates  [12]. 
Clearly,  additional  analysis  of  the  effects  of  curvature  on 
composite  materials  Is  required. 


SOLID  PANEL 


NATURAL  FREQUENCY 


VI  I .  CONCLUSIONS 


In  addition  to  the  observations  and  conclusions  In 
previous  chapters,  the  fol lowing  additional  conclusions  are 
i ncl uded: 

1.  The  STAGS  C-1  computer  code,  assuming  small 
displacements  and  linear  stress-strain  relationships,  proved 
excellent  at  predicting  the  mode  shapes  o-f  the  composite 
panels  tested  here. 

2.  Similiarly,  the  STAGS  C-1  code  proved  capable  of 
predicting  natural  frequencies  within  limits  of  engineering 
accuracy. 

3.  For  the  c 1 amped-c 1 amped  boundary  conditions,  the  QUAF 
410  plate  element  proved  as  effective  as  the  higher  order 
QUAF  411  el ement . 

4.  For  the  smallest  cutout,  the  mode  shapes  change  very 
little  in  comparison  with  the  solid  panel. 

5.  Natural  frequencies  tend  to  decrease  with  cutout  size. 


6.  Mode  shapes  for  the  larger  cutouts  are  significantly 
altered.  The  Inclusion  of  mode  shapes  not  found  In  the 
solid  panel  as  well  as  symmetry  shifts  In  fundamental  modes 
were  observed  for  Increasing  cutout  size. 


RECOMMENDATIONS 


1.  Using  both  finite  element  analysis  and  holographic 
interferometry,  the  range  at  this  study  should  be  extended 
to  higher  frequencies  and  mode  shapes.  Similiarly,  the 
study  on  nonsymmetric  laminates  should  be  included. 

Verifying  the  accuracy  of  the  finite  element  code  to  model 
the  coupling  of  bending  and  extensional  stresses  should 
extend  confidence  in  the  code  for  these  type  of  laminates. 

2.  Interchanging  the  order  of  the  45  degree  lamina  should 
prove  instructional  on  their  effect  on  the  left  tilt  seen  in 
the  displacement  patterns. 

3.  The  analytical  techniques  used  in  this  study  should  be 
extended  to  other  boundary  conditions. 
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SAMPLE  INPOT  FILE  2X4  INCH 

2,0, 1 ,0,0, 0,0, 0,0 

1,0,0 
1,0, 1,0 


1,0,7000,1 
5, 0,0,0 

25,25 

1,0 

1 8. 84E06, .021813, .9099E06, 

1,1,8 

1 . . 005.0. .0 

1.. 005,-45.,0 

1. . 005.45. .0 

1.. 005.90..0 

1.. 005.90..0 

1.. 005.45..0 

1 . . 005,-45. ,0 

1 . . 005.0.  ,0 
5,0 


0.0,12.0,-28.648,28.648,12. 
XI , X2 , Y1 , Y2 (DEGREES ) , RADIUS 
1,0 

410,0,0,1,0,0,0 

11,15,11,15 
2, 2, 2, 2 
0 
1 


CUTOUT  410  BLMT  $A1 

$B1  SMALL  VIBRATIONS, SAVE 
MODEL  DATA, GRID  TEST 
$B2  ONE  SHELL 
$B3  SINGLE  MATERIAL, ONE 
WALL  TYPE 

SCI 

$D2  7000  CPU  SECONDS 
$D3  EIGENVALUES, SHIFT, UPPER 
LOWER  BOUNDS 
$F1  ROWS, COLUMNS 
$11 

055,. 0216216, 1.468E06, 15. 2  $12 
MATERIAL  PROPERTIES 
$K1  GENERAL  LAYERED  WALL- 
EIGHT  LAYERS 

$K2  PLY  THICKNESS  -  FIBER 
ANGLE 

$K2 

$K2 

$K2 

$K2 

$K2 

$K2 

$K2 

$N1  CYLINDER  LOCAL 

COORDINATES  COINCIDE 
WITH  GLOBAL 

$N2A 

$M5  NO  IMPERFECTIONS 
$N1  410  ELEMENT  25  X  25 
GRID 

$N4  CUTOUT  DEFINITION 
$P1  CLAMPED  ALL  EDGES 


$R1  PRINT  OUT  DISPLACEMENTS 


I  4 


I  J 

••  «•  V 
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APPENDIX  A 


SAMPLE  INPUT  FILE  SOLID  PANEL 

2,0,1 ,0,0, 0,0, 0,0 

1r0,0 
1,0, 1,0 

1 

1,0,7000, 1 
5, 0,0,0 

25,25 

1,0 

1  8. 84E06, .021813, .9099E06,. 

1,1,8 

1 . . 005.0. .0 

1 .. 005,-45. ,0 

1.. 005.45..0 

1. . 005.90. .0 

1.. 005.90..0 

1.. 005.45..0 

1.. 005,-45.,0 

1 . .  005 . 0 . . 0 
5,0 


0.0,12.0,-28.648,28.648,12. 
X1,X2,Y1,Y2( DEGREES), RADIUS 
1,0 

410,0,0,0,0,0,0 

2, 2, 2,2 
0 
1 


410  ELNT  $A1 

$B1  SMALL  VIBRATIONS, SAVE 
MODEL  DATA, GRID  TEST 
$B2  ONE  SHELL 
$B3  SINGLE  MATERIAL, ONE 
WALL  TYPE 

$C1 

$D2  7000  CPU  SECONDS 
$D3  EIGENVALUES, SHIFT, UPPER 
LOWER  BOUNDS 
$F1  ROWS, COLUMNS 
$11 

,. 0216216, 1.468E06, 15. 2  $12 
MATERIAL  PROPERTIES 
$K1  GENERAL  LAYERED  WALL- 
EIGHT  LAYERS 

$K2  PLY  THICKNESS  -  FIBER 
ANGLE 

$K2 

$K2 

$K2 

$K2 

$K2 

$K2 

$K2 

$M1  CYLINDER  LOCAL 

COORDINATES  COINCIDE 
WITH  GLOBAL 

$N2A 

$N5  NO  IMPERFECTIONS 
$N1  410  ELEMENT  25  X  25 
GRID 

$P1  CLAMPED  ALL  EDGES 

$Q1 

$R1  PRINT  OUT  DISPLACEMENTS 
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